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A density- functional formalism for superconductivity and magnetism is presented. The resulting 
relations unify previously derived Kohn-Sham equations for superconductors and for non-collinear 
magnetism. The formalism, which discriminates Cooper pair singlets from triplets, is applied to two 
quantum liquids coupled by tunneling through a barrier. An exact expression is derived, relating 
the eigenstates and eigenvalues of the Kohn-Sham equations, unperturbed by tunneling, on one side 
of the barrier to the proximity-induced ordering potential on the other. 
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I. INTRODUCTION 

Condensed-matter phenomena owe much of their variety to the multifarious properties of inhomogeneous electron 
liquids. Long-range order, gaps, screening, enhanced correlations, anomalies, resonances — a multiplicity of effects 
results from the diversity of chemical compositions and crystal structures. Not surprisingly, density-functional the- 
ory (DFT), a programme dedicated to describing the effect of atomic-scale inhomogeneity on the electron liquid, 
has acquired mounting prominence. ^'^ Its success in the description of properties associated with the microscopic 
inhomogeneities due to lattice potentials provides the foundation upon which ah initio procedures are built. To 
describe ordered states, specialized formalisms comprising other order parameters along with the density have been 
proposed. "^"^ 

Macroscopic inhomogeneities have received less attention from the DFT community. In particular a physical 
boundary separating distinct order parameters is a nesting ground for concepts and applications, of which the discovery 
of giant magnetoresistances is the most persuasive example.*'^ More recently, advances in fabrication techniques and 
experimental probes (notably X-ray and neutron reflectometry) have revealed that the behaviour at such boundaries 
can be quite unexpected. To give one example, novel magnetic order has been observed in the superconducting side 
of a ferromagnet-superconductor interface [Ref. 10; see also Ref. 11 and references therein]. Yet one finds in the 
literature no DFT broadly applicable to junctions. 

To fill this void, we present here a density-functional theory of superconducting and magnetic materials. A gener- 
alization of previously published theories for the individual orderings,'^"^ the new formalism is designed to describe 
long-range order with coexisting charge, magnetic and superconducting order parameters. Competition or co-existence 
of different forms of order is known to occur in bulk high-temperature,^^ heavy fermion,^^'^* and organic^^ super- 
conductors as well as manganites.^^ Here we will focus on proximity effects in junctions between differently ordered 
quantum liquids. 

While studies of superconductivity in a specific material can be restricted to Cooper pairs of given spin, either 
singlets^'^^ or triplets, the more general setups that we target require parallel treatment of the singlet pair, the 
triplet pair, and the magnetization densities. A triplet superconductor coupled to a non-collinear antiferromagnet is, 
e.g., well within the scope of our formalism. 

As a general application, we consider two quantum liquids with different order parameters separated by a thin barrier 
that allows tunneling. Under these circunstances, we demonstrate that the Kohn-Sham (KS) equations yielding the 
ground-state energies and densities on one side of the barrier can be decoupled from the analogous equations for the 
energies and densities on the opposite side. In each decoupled KS Hamiltonian, an effective potential obtained from 
the solution of the unperturbed KS equations, i. e., the KS equations in the absence of tunneling, represents the 
opposite side. This potential gives mathematical substance, within DFT, to the proximity effect. In a normal metal- 
superconductor junction, the decoupled KS Hamiltonian adds correlation to the Bogolubov-de Gennes equations. ^^'^^ 
In a normal metal-antiferromagnet junction, it generates analogous equations describing the staggered proximity field 
induced on the nonmagnetic side. 

Our presentation starts out with a cursory review of DFT for superconductivity and for magnetism. Section III 
presents the formalism for coexisting order parameters. The resulting KS equations are applied in Section IV to a 
barrier separating two quantum liquids, and the effective proximity potential is derived. In Section V the formalism 
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is generalised to magnetic interfaces. Finally, Section VI lists our conclusions. 

II. DENSITY-FUNCTIONAL THEORY FOR ORDERED QUANTUM LIQUIDS 

The following two subsections briefly recapitulate the DFT approach to ordered quantum liquids. They highlight 
those aspects of the theory that will prove particularly important in our formulation and others that were brought to 
view in recent publications, posterior to the original references. 

A. Density-functional theory for superconductors 

Superconductors still pose a challenge to electronic-structure theorists. Progress in that area had a late start. 
Long after band-structure calculations based on DFT had provided valuable information about such normal-state 
properties as Fermi surface geometries, single-particle spectra, and electronic densities of states, well after DFT had 
yielded such lattice properties related to superconductivity as phonon dispersion relations, the superconducting state 
still lay outside the realm of ab initio electronic-structure calculations. 

Model Hamiltonians had provided much of what was known about superconductivity. The reduced BCS 
Hamiltonian,^^ the Hubbard Hamiltonian^^ and its variations, and the Bogohibov-de-Gennes mean-field equations"'^^'^'^ 
are examples, none of which seemed adaptable to a density- functional formulation. In view of its reliance on empir- 
ical information on phonon spectra and Coulomb matrix elements, not even the detailed microscopic description of 
strong-coupling superconductivity in Eliashberg's theory could be seamed to ab initio DFT.^^ 

To circumvent such difficulties, an alternative approach was proposed.^ Instead of functionals of the density, one 
now studied functionals of two variables: the normal density n(r) = (^'j(r)^'|(r)) and the superconducting order 

parameter x(r,r') = (^|(r)*I'|(r')). In the same way that p is coupled to an electric potential, or the magnetization 
density of spin-DFT (SDFT) is coupled to a magnetic field, the anomalous density x was coupled to a pair potential. 
In the same way that the KS equations generalize the Hartree mean-field equations, extended KS equations were 
derived that generalize the mean-field Bogolubov-de Gennes equations. 

Two are the potentials in these KS Bogolubov-de Gennes equations: an eff'ective electric potential 

Vs{r) = Vext{r) + VH{r) + Va;c{r), (1) 

and an effective pair potential 

A(r, r') = Ae,t(r, r') + (r, r') + A,e(r, r'). (2) 

If the (microscopic) inliomogencity is due exclusively to the lattice potential v^xt, only an infinitesimal external 
potential Aext is needed, to break the gauge symmetry that would otherwise annul the anomalous density in the 
self-consistent cycle. (As pointed out in Ref. 6, however, and further discussed in Section IV A, macroscopic inho- 
mogeneities may generate non-infinitesimal external anomalous potentials.) While the Hartree potential Vh tends 
to make the charge distribution uniform, the anomalous Hartree potential Ajf, the interaction of the anomalous 
density with its own pair potential, tends to enhance the superconducting order parameter. As usual, the exchange- 
correlation potential Wa;c(r) = SExc[n, x]/Sn{r) is the derivative of the universal xc functional of DFT with respect to 
the normal density. Similarly, A2;c(r,r') = SExc[n,x]/^x{'<^T'<^') represents the exchange-correlation correction to the 
mean-field approximation. Among the generalizations of this formalism, we mention one that will assist our analysis: 
the extension to triplet superconductors.^*^ 

The exchange-correlation functional Exc[n,x] is, of course, imknown, and the DFT programme calls for first- 
principles approximations. That at least in the context of plionon-mediated superconductivity this programme can be 
followed to its end was demonstrated by Gross and collaborators, ^^"^"^ who constructed a functional with no adjustable 
parameters and applied it to a number of materials. We note in passing that their breakthrough has led to the first 
truly microscopic theory of conventional superconductivity. 

B. Density-functional for magnetic systems 

The fundamental variables in the standard formulation of DFT, collinear spin-DFT (SDFT), are the spin-resolved 
densities ni-(r) and nj,(r). From these, the charge density n(r) = ni-(r) -|- n|(r) and the magnetization mz{r) = 

/i(n|(r) — ^^(r)) are promptly recovered (/io is the Bohr magneton). We concentrate our discussion on two features of 
magnetic systems that are foreign to SDFT in the local spin-density approximation: non-collinearity and nonlocality. 
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To describe non-coUinear magnetic structures, one may substitute the magnetization vector ■m(r) for its z- 
component. Elegant extensions of non-coUinear SDFT and succesful implementations have been cast in this forge. ^-"^^^^ 
However, while local approximations suffice to describe the nearly uniform average local magnetization in a fcrromag- 
net, such approximations cannot be expected to fully reproduce the strong nonlocal correlations associated with spin 
waves in a non-collinear antiferromagnet.^^ 

Faced with this difficulty, one might construct nonlocal, e.g., orbital-dependent, functionals of the local variables 
n(r) and m(r).^^ A simpler alternative is suggested by the procedure that, starting with spin-independent DFT, 
constructed SDFT and DFT for superconductors: a density sensitive to the ground-state correlations characteristic 
of the phenomenon under study is added to the list of fundamental variables. To the set {ni-(r), nj^(r)} of SDFT 
variables, non-collinear antiferromagnetism thus adds the staggered density^ '^"^ '^^ 

Ps{v,v') = {i>]{v)i>^{v')). (3) 

The formal similarity between ps and the anomalous density in DFT for superconductors is not accidental: already 
in 1958 it was recognized that the restricted particle- hole transformation^^'^° 

- *l(r) 

converts the BCS procedure into a mean-field theory of antiferromagnetism. 

Central to the development of staggered DFT is a coupling between the staggered density and a staggered potential 
^(r, r'), a nonlocal generalization of the magnetic field B(r). In the resulting KS equations, the staggered density is 
coupled to the effective staggered potential 

Ss{t, r') = Sext (r, r') + ^^(r, r') + 5,c(r, r'). (4) 

Exploring the analogy with the external pair potential in DFT for superconductivity, Refs. 7 and 37 conjectured that, 
more than a mathematical artifact, the first term on the right-hand side could be interpreted as a proximity effect, a 
potential induced near a non-collinear antiferromagnet. 

Two features distinguish staggered DFT from collinear SDFT. (i) The diagonal element (r' = r) of the staggered 
density determines the x- and y-components of the magnetization: mx{v) ~ /io[/0.s(r, r) + /7*(r, r)] and ray{v) = 
i/io[Ps(r,r) — Ps{v,y)]. This restriction yields a formalism equivalent to local non-collinear SDFT.'^^""^^ Even under 
this restriction, even for perfectly collinear states, non-collinear DFT is more powerful than standard (i. e., collinear) 
SDFT. In the latter, the specification of the quantization direction breaks rotational symmetry. The staggered 
density restores that symmetry: collinear SDFT is blind to, e.g., a magnetization along the x axis, since = = 0. 
Staggered SDFT, by contrast, extracts the x and y components of the magnetization from pg. 

(ii) The nonlocal dependence Ps = Ps{y,y') enhances the superiority of staggered DFT over SDFT, and makes it 
more powerful than local non-collinear SDFT. To see this, it is sufficient to recall that the electron liquid possesses an 
instability against the formation of spin-density waves (SDW).^^~^^ The driving force of the Overhauser instability is 
the staggered Hartree interaction 

UM = -jd'rjdS'\B^^^. (5) 

This interaction, which formally appears as a Hartree term involving the staggered density, arises from evaluating 
the exchange diagram with two-component spinors. It tends to push the energy of the SDW state below that of the 
paramagnetic state. '^'^'^ 

Standard formulations of SDFT miss this energy reduction altogether; sophisticated approximations to the 

exchange-correlation functional are needed to account for it.^^ By contrast, the right-hand side of Eq. (5) is trivially 
incorporated in the definition of the staggered DFT exchange-correlation energy functional; the epithet "staggered 
Hartree energy" emphasizes its formal similarity to the usual Hartree term while the subscript x is a reminder of 
its physical origin in the exchange diagram.^'^^'^^ Its contribution S'_f/(r,r') = 5Ux/ 5ps{jc,v') to the effective Hartree 
potential 5*8 (r, r') defines the second term on the right-hand side of Eq. (4). Staggered DFT has been shown to recover 
both the Overhauser instability in the exchange-only approximation and its suppression by correlation.^ 



III. DENSITY-FUNCTIONAL FORMALISM FOR COMPETING ORDER PARAMETERS 



We now turn to a formulation that extends the results recapitulated in Section II to more complex ground states, 
with coexisting or competing order parameters. Our derivation following closely Ref. 18 for the superconducting order 
parameter and Ref. 7 for the non-collinear order parameter, we list but the key equations. 
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A. Many-body Hamiltonian 

We consider an interacting Hamiltonian that includes external potentials coupled to all the physical observables of 
interest: 

H = f + U + [ d^r flair) va{r) 

a- ^ 

+ / (i'x(fy{ps{x,y)S{x,y)-x{x,v)^.*{x,y)- ^ x^{x,y)l^l,{x,y) +11. c). (6) 

•' m=-l 

Here, the density operators na{r) = '^)^{r)'^cT{i'), ps{r,r') = X(^)^') = [V'(^)V'(^')]s i^^^ singlet com- 

bination of the two field operators), and x„(r,r') = [ip{r)il;{r')]m (the m = — f,0 and 1 triplet combinations of the 
field operators) are coupled to the external potentials v^r {r) = v (r) + afisB^ (r), S{x,y) (the staggered potential^) 
A(x, y) = A(y, x) (the singlet pairing potential^^), and A^ix, y) = — A„i(y, x) (the rn component of the triplet pair- 
ing potential -'^®), respectively. As usual, the kinetic energy operator is T = X)<t=t i /'^'^^'^i (^) (— ft^V^/2m) ip^ (r) 
where m is electronic mass, and as in Refs. 6,f8, we assume the following form for the electron-electron interaction: 

^■=\ X] j d^y\d^xid^X2d^y2'ipl^{vi)i)l,^{xi)U {yi,(3i]Xi,ai]X2 (7) 

/9i,Qi,a2,/32 

where the electron-electron interaction comprises the Coulomb repulsion and a nonlocal spin-dependent term, needed 
to represent, e.g., the phonon- mediated attraction in conventional superconductors or a spin- fluctuation mediated 
interaction in unconventional ones: 

U {yi,(3i;xi,ai;x2,a2;y2,/32) = i 7S{xi - X2)S (yi - y2) SauaJ^u^^ + W {yi, (3i;xi,ai;x2,a2;y2, P2) ■ (8) 

\xi - yi\ 

Since the following arguments make only implicit reference to W, to be concise we have chosen a non-retarded potential 
to illustrate Eq. (8). 

Lest the reader be puzzled by the asymmetry in our treatment of the pair density, which is resolved into a triplet 
and a singlet components, while the staggered density is not, we note that the restricted particle-hole transformation 
V't(r) i'li'i'), ipi{r) I— > 4'\{'>') would turn spins into isospins. Under this transformation, the Cooper pairs would 
comprise an isospin dublet coupled to an anomalous magnetic-field like potential, while the staggered density would 
decompose into an isospin singlet and an isospin triplet. Although the two approaches are mathematically equivalent, 
we find the language of spins more attractive for the present purposes than that of isospins. 

B. Kohn-Sham equations 

Given U, the ground-state energy is a functional of the densities 71^ {a =t, J,), Ps, X> and Xm ("^ = ~1)0, 1), which 
can be written as 

E[n^,ni,Ps,X, Xm=0 > Xm=+1 > Xm= - 1] = ^H^f [«-T > > > X ' Xm=0 > Xm=-|- 1 > Xm= - 1 ] + 

J2 [ d^rn„{r)v„{r)+ j d^xd^y{ps{x,y)S{x,y) - x{x,y)A*{x,y) - ^ Xm{x, y)^*m{x, v) + cc), (9) 

<T '' m=-l 

where the potential energy in the various external potentials has been written explicitly, and the kinetic and interaction 
energy are combined into the Hohenberg-Kohn internal-energy functional Fhk '■ 

Fhk = (f) + {U) =T, + Umf + E^^. (10) 

In the last equation we defined the exchange-correlation functional E^c = (T) — T, + (U) — Umf in terms of the 
kinetic energy of a noninteracting system with densities ni,ni, pa,x,Xm=0'Xm=+ijXm=-i' ^^'^ the mean-field 
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approximation Umf to the full interaction energy {U). Formally, this mean-field approximation is given by 



/3i,a;i,a;2,/32 

f)i,ai,a2,02 

j d^yid^xid^X2d^y2{'ipl^{yi)i>a2{x2))U {yi,pi;xi,ai;x2,a2-,y2,p2) {■tl>i,{^^^ (11) 

/3i,ai,Q2,/32 

Some contributions to Umf have appeared before in more restricted formulations, such as the conventional Hartree 
term 



E, = Ud^x[d^y^p^^ 
2 i J |x - y 



the anomalous Hartree term 



and the staggered Hartree term 



E'„=J d'x j 



|x-y| 



E!, = -ld'xl d\ ^^^^^ = UAPs\. 



(12) 



(13) 



(14) 



We stress that while E'^^ is repulsive, and thus detrimental to superconductivity,^'' t/^iPs] is attractive, and favours 
formation of spin-density wavcs/'"'^''^** 

The KS procedure, applied to this total-energy functional, leads to single-particle equations for the four-component 
spinor 



(15) 



with particle (hole) wavefunctions Ua{r) [va{r)] for each electron spin, a =t, i- The symbol T indicates transposition, 
while the superscript £ labels the self-consistent solutions of the KS equations, the first of which is 



w 



(16) 



where w^^^ denotes the set of the seven potentials {vi^\S<^^\A^'\A^r^^} {a =t, i, m = -1, 0, 1). The KS Hamiltonian 



IS 



2A 



V a[,")* - AW* 2AL1* 



.5(^)5 



(17) 



with the shorthands h for —h^V"^ /2m, and S^-'^\ A^"^ and A^^ for the integral operators associated with the non-local 
potentials (r,r'), A(^) (r,r') and A^^' (r,r') by 



5W: /(a;) ^ j d^yS^^\x,y)f{y), 

A(^) : fix) ^ J d^yA^'\x,y)f{y) and 

AW: /(x) ^ y d32/AW(a;,y)/(y). 



(18a) 
(18b) 
(18c) 
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The effective single-body potentials are defined by 

Vas = W<T + — F h -z , (19a) 

Oitfj Oit(j 

AW=A+^^ + ^, and (19c) 

ox ox 

"Am "Am 

The derivation of Eq. (17) from Eq. (9) is analogous to that of the simpler KS Hamiltonians in Refs. 6 and 7. First, 
we note that minimization of the interacting Hamiltonian of Eq. (6) with respect to results in the Euler equation 

+ Va{Y)=Q, (20) 



6na{r) (5na(r) 6na{v) 

and similar equations for the minimization with respect to the other densities. 

In a second step, we construct a non-interacting many-body system with Hamiltonian Hg obtained by substracting 
IJ from the interacting Hamiltonian of Eq. (6) and replacing the external potentials ii^, S, A and with u^-.s) S'^^\ 
A^'^) and Am"*, respectively. These self-consistent potentials are chosen such that the non-interacting system has the 
same densities n-f, n^, ps, X and Xm as given by the original Hamiltonian, H. Minimization of Hg with respect to 
leads to the Euler equation 

+t;.,«(r)=0. (21) 



The two Euler equations can be combined to yield Eq. (19a). In the same manner, minimization with respect to 

the other densities yields Eqs. (19b,19c,19d). To solve the effective problem, we need to diagonalise Hg- We achieve 
this by a a linear (Bogolubov) transformation from the electron creation and annihilation operators, V'^C'') and V'o-Cr), 
to new fermionic creation and annihilation operators , 7; : 

^,{v)=Y,[<{^)li+v'A^ri}\. (22) 

This transformation renders the effective single-particle Hamiltonian in diagonal form: Hs = constant -|- ^lljll- 
The single-particle energies ei are the eigenvalues of the KS Hamiltonian (17), while the coefficients wj^ and v^^ of the 
Bogolubov transformation arc the components of the spinors 

In terms of these coefficients, the ground-state densities are given by 

(23a) 
(23b) 

and (23c) 
(23d) 

The ground-state energy is then obtained by computing the functional £[71-^ ,ni, ps, X) Xm=05 Xm=-i-ii Xm=-i]- 

The KS Eqs. (16-19d) determine the ground-state energy and densities for external potentials ranging from micro- 
scopically to the macroscopically inhomogeneous, as well as from periodic to disordered. Our interest in the proximity 
effect focuses the following discussion on the inhomogeneities associated with a one-dimensional discontinuity sepa- 
rating two semi-infinite regions. 



ria- (r) 


= E 


H{rf+ E l<(^)r (^=T.i) 




i:ei<0 




Ps{x,y) 


= E 


u\*{x)ul{y)+ v'^{x)vl*{y), 








x{x,y) 


= E 


[v'*{x)u'{y)l+ J2 [uHx)v'*{y)U 




i:se<a 


e:ee>0 


Xm{x,y) 


= E 


[v'*{x)u\y)U+ Y [u'ixy*iy)]m. 






t.6e>0 
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Figure 1: The physical system under consideration (see text): the shaded area represents the physically inaccessible region 
across which tunneling can take place, while the dashed lines indicate the extent of the L and R regions on which the tunneling 
amplitude t (r, r') is important. 



IV. SPECIALIZATION TO INTERFACES 



We consider two semi-infinite spatial regions, labelled L (left) and R (right), separated by a potential barrier that 
is wide in comparison with the range of the interaction U in either L or R. Though high, the barrier is finite and 
hence allows electronic tunneling between L and R. As Fig. 1 indicates, we denote the position in L (R) by r (r'). 
To save space, we will indicate the densities and potentials on the R semi-space by a prime. Thus, e.g., v{r £ L) 
[v{r e R)] will be denoted v [v'). 

Our model Hamiltonian reads 



H = + Ht + Hf. 



(24) 



Here, Hl [Hr) are of the form H + V, where H is defined in Eq. (6), and V is an infinite barrier centered at the 
interface between L (R) and the tunneling region, which enforces orthogonality between the eigenfunctions of Hl and 
Hr.'^^ The phenomenological tunneling Hamiltonian 



3„/ 



„ Jl Jr 



ipl (r) t{r,r') i}„{r') + H.c. 



(25) 



provides transport across the barrier. 

We complete the phenomenlogical description of the interface by assuming that the only effect of the tunneling 
term Ht in the Hamiltonian (24) is to make the KS equations become 



H, 



(s)t 




(26) 



where the tunnehng KS Hamiltonian is defined by 



















i 
















-i* 







Vo 








-i* 


1 



(27) 



The integral operator i is given by t : f [r) ^ J Sr' t (r, r') f (r'). 

Our phenomenological description of the barrier contrasts with the microscopic treatment of the semi-infinite re- 
gions. Instead of deriving and solving Kohn-Sham equations for the shaded region in Fig. 1, we prefer to follow 
Prange's prescription'*'*, even if we have to rely on phenomenological considerations to define the tunneling matrix 
t (r, r'). This simplification focuses our analysis on the interaction between the L and R regions. Under appropriate 
circumstances, as shown below, it leads to an independent set of Kohn-Sham equations for each region, a result com- 
bining mathematical convenience with physical appeal. Under such circumstances, to determine the densities in one 
of the regions, we can rely on first-principles calculations in which the external potential captures phenomenologically 
the influence of the material across the barrier. 
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A. Decoupled Hamiltonians 

The KS equations (26) describe the electronic state of L and R. Let us now derive a set of KS equations that 
describe, formally, only the states of R, but take into account the influence on such states of the charge, spin and 
superconducting order that are potentially present in L. 

For r G L, since the eigenfunctions of the KS Hamiltonian with Ht = form a complete set, we can write 

*L = E(^l^*^)*i' (28) 
e 

where the tildes denote the absence of tunneling, and the scalar product ^^i,^^^ comprises integration over the 

spatial variable and summation over spinor indices. 
Insertion into the second row of Eq. (26) then yields 



(29) 



To eliminate $^ from this equation, we first write the eigenvalue equation defining i. e., the first row of Eq. (26) 



with H):' = 0, 



w 



$1 = ^1*1, 



and multiply it on the left by . Hermitean conjugation then yields 



w 



is) 



Next, we multiply the first row of Eq. (26) by on the left: 



(30) 



(31) 



(32) 



We then subtract Eq. (31) from this result, solve for l>^^$^, and integrate over the spatial variable to find that 

(^l^^^) = ^^K^^^'^«)• (33) 

We can now substitute the right-hand side for the scalar product on the left-hand side of Eq. (29) and exploit the 
completeness of the wavefunctions (Z = 1, 2, . . .) to define the effective external Hamiltonian 



(34) 



so that Eq. (29) defines a KS eigenvalue equation for r' G R: 



w 



(«)' 



(35) 



More explicitly, given that, if (Mf, U|, w-f, {)|)^ is an eigenvector, with eigenvalue e, then (u-f, Uf, uj^)^ is also an 
eigenvector, with eigenvalue — e, we can cast the effective external potential in the form 

/ vi^'^ie.) 5W(£.) -2A(^)(e,) -(A^^) + A^^^) (e,) \ 



t 



-2AL1 (e. 



S^'^He.) i?{e.) (A(^)-A(^))(e.) 

2A(^)*(-e.) (A[/)*+A(^)*)(-e.) -i)^ (-£.) -S(^> {-e.) 

[{Ai'> -A(^>){-s,) 2ALt(-£.) -5(«)^(-£.) -i>['H-£.) J 



(36) 



Figure 2: The external non-local "potential" of Eq. (37a) induced in the R region describes new coupling between points x' and 
y' that takes into account processes in which an electron would tunnel across te barrier, propagate inside the differently-ordered 
electron liquid in the L region, and then tunnel back into R (see text). 



where the integral operators {su), S {e„), A {e„), and (ej,) (m = —1, 0, 1) have the following kernels: 

v,{s,;x',y') = [ d^xt{x',x) V / t{y,y'), (37a) 

Jl i £iy-£Ld 

S{6,;x',y') = [ d^xtix',x) V / <T(^)<i(^) tiy,y'), (37b) 

Jl e ~ ^^'i' 

A{s^;x',y') = [ d'xt{x\x) V / rf^y M^^^Ml {-l)t*{y,y'), and (37c) 
Jl ^ Jl — £l,i 

Am{e.;x',y') = [ d'xt{x',x) ^ / MM^|Mk (-l)t*{y,y'). (37d) 

The effective external potential H^xt {su) expresses mathematically the proximity effect. The non-local, effective 
potential {£v',x' ,y') in Eq. (37a) describes the virtual transition from point x' to y' depicted in Fig. 2, with 
intermediate tunneling to point x and propagation to point y. The right-hand sides of Eqs. (37b-37d) have analogous 
interpretations, the propagation in L now involving a spin flip [Eq. (37b)] or spin-conserving scattering into a hole 
state [Eqs. (37c) and (37d)], respectively. 



B. General aspects of the decoupled Hamiltonian 

Equations (37a-37d) decompose a pair of quantum liquids, differently ordered and coupled by tunneling, into two 
uncoupled liquids subject to effective, energy-dependent external fields. The right-hand side of each equality depends 
only on not on or on In words, the effective fields induced by the proximity effect on one half-space 
depend only on the solutions of the KS equations on the opposite side, unperturbed by tunneling. 

In general, the effective fields also depend on the energies e^. Whether it is easier to diagonalizc energy independent 
coupled Hamiltonians or energy-dependent uncoupled ones then depends on specific aspects of the problem under 
study. If a low barrier separates degenerate states, the coupled-Hamiltonian formalism is natural and mathematically 
convenient. As pointed out in Ref. 6, however, if the energy spectrum e on one of the sides is bounded by a gap, while 
that on the other side is not, then the low KS eigenvalues can be neglected in the denominators of the sums on the 
right-hand sides of Eqs. (37a-37d), and the potentials become energy independent. In this case, the treatment of the 
decoupled Hamiltonians is clearly more efficient than the solution of the coupled KS equations. Prom a conceptual 
viewpoint, moreover, the emergence of proximity-induced effective potentials that give rise to superconductivity or 
antiferromagnetism is appealing. 

In the analysis of infinite systems, at the formal level, DFT starts out by introducing a ficticious external potential, 
set equal to zero after the KS equations are derived. Equation. (6) is an example. To solve the KS equations, at the 
operational level, one likewise seeds the self-consistent loop with an artificial external potential. In the vicinity of an 



10 



ordered system, as Eqs. (37a-37d) show, this mathematical expedient is neither formally nor operationally necessary. 
One can therefore interpret the spontaneous ordering of an infinite system as a residual consequence of a proximity 
effect. We see that this interpretation, first offered in connection with superconductors in Ref. 6, describes equally 
well magnetic, superconducting and coexisting states. 

In spite of the complex form on the right-hand side of Eq. (36), our attention to spin in the definition of the 
Cooper pair densities guarantees that the KS equations couple only the singlet (triplet) density to the singlet (triplet) 
potential. As the analyses in subsections IV C and IV D will show, if a conventional superconductor is coupled to a 
non-collinear antiferromagnet by an interface through which electrons can tunnel, the resulting proximity potential 
will induce the formation of singlet Cooper pairs in the antiferromagnetic material; likewise, on the superconducting 
side, the induced staggered field scatters the Cooper pairs into other singlet states. 



C. Special case: superconductivity 



This subsection and the next examine Eqs. (37a-37d) in two particularly simple situations. Here, we detail the first 
arrangement, which couples an unspecified material (R) to a singlet superconductor (L). We take the latter to be 
nonmagnetic, so that Vs,i — = Ss = 0. 

To determine the unperturbed KS eigenstates in L, which contribute to the right-hand sides of Eqs. (37a-37d), we 
set = 0. A Bogolubov-Valatin transformation then diagonalizes the matrix on the right-hand side of the resulting 
Eq. (17). The eigenspinors have the form 



( <T(r) \ 




(5„,2 (a = 1,2) 



(38) 



The amplitudes u^(r) and vz{y) are solutions of the Bogolubov-de Gennes equations [i.e., Eq. (7) in Ref. 6]. The 
invariance of the L Hamiltonian under the inversion z ^ —z makes the KS eigenvalues independent of the spin index 
a, i.e., ejj = e^. The eigenvalues moreover change sign under the transformation (ui, Uj, Uj)"^ hh. (u|, W|, {t|, -uj^)^. 
These two symmetries considered, and the eigenfunctions on the right-hand side of Eq. (38) substituted on the 
right-hand sides of Eqs. (37a-37d), we find that 



Jl , Jl en 



A{e^;x',y') = / d^xtix',x)y] / d^y 
Jl jl Jl 

S{£„;x',y') = Am{£u;x',y') = (m 



ue {x) v} (y) 



El — 
-1,0,1). 



t* {y, V') , and 



(39a) 

(39b) 
(39c) 



Thus the superconductivity in L induces only a normal and a singlet pairing potential in R. When i? is a semi- 
infinite, normal metal, we recover the SC-DFT result.^ Equation (39a) is the nonlocal normal potential alluded to, 
but not given explicitly, in footnote 10 of Ref. 6. 



D. Special case: magnetism 



As a second particular case, we consider a non-superconducting material in L, i.e., A/, = 0. The eigenspinors of 
the KS equations now have the form 



$r(r) 



5™(r) 









(t =p, h), 



(40) 



where the label t = p and t = h designates particle- and hole- like KS quasiparticles, respectively. The amplitudes 
^m,a{^) are solutions of an eigenvalue problem analogous to the Bogolubov-de Gennes equations [Eq. (15) in Ref. 7]. 
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For a given m, the corresponding eigenvalue can either be positive or negative, and the external potentials are 

(e.; x\ y') = / d\ t xjj^l d'y ^--^ ('"X^ t (y, y') (41) 

S{e.;x',y') = [ Sxt{x',x) ^ / d'y ^I^^:l^^±^kMt^y,y') (42) 

Jl „ J L — Em 

A(e,;a;',y') = A™(e,;a^',y') =0 (m = -1,0,1). (43) 

This set is analogous to Eqs. (39a-39c) and gives substance to the image^ of the external contribution to the staggered 
field in DFT as a proximity eflfect, an interpretation that places additional emphasis on the analogy between DFT for 
superconductors and for spin-density waves. 



V. GENERALIZATION: MAGNETIC INTERFACES 



The tunneling Hamiltonian (27) defines a inert barrier, one that merely allows charge transport between regions L 
and R. By contrast, this section examines a magnetically active interface, in which spin-orbit coupling or spin-flip 
scattering from impurities add off-diagonal elements to the tunneling matrix. In special situations, a trivial rotation 
of the quantization axis may be sufficient to diagonalize that matrix over the entire barrier; in such instances, one 
expects the barrier to induce spin polarization in both L and R. 

More generally, however, the spatial dependence of the matrix elements will bar global diagonalization, so that 
instead of simply redefining the quantization axis, the barrier will turn into an inhomogeneous source of spin fiips. 
While the system, constituted by the L and R regions and the interface between them, is subject to global conservation 
laws, the spin of the electrons in regions R and L is no longer conserved. This section explores the consequences of 
that rupture. We shall see that while the potentials v^, {a =|, |) and S are only quantitatively affected, the magnetic 
barrier breaks the independence between the singlet and the triplet Cooper pairs. A magnetic interface separating a 
singlet superconductor from a normal metal, for instance, will induce the formation of triplet pairs in the latter. 

More specifically, instead of Eq. (27) we consider the following tunneling KS Hamiltonian 



Tmagn 



/in hi \ 

tiT hi \ \ 





(44) 



V -ti^ -ti^ J 

and follow the analysis leading from Eq. (27) to Eq. (34), which now takes the form 

Hext,magn {^u) = X] ^rliagn^L ^ _ ~l ^"^LHrliagn- (^5) 



From Eqs. (15) and (44), the computation of the right-hand side of Eq. (45) is straightforward, which brings Eq. (45) 
to the form of Eq. (36), with the kernels 

Vr{e.;x' ,y') = [ d^x J^trA^^x) [ d% Y^ O^lAv) t^,^^y^y')^ (r=T,i) (46a) 

Jl Jl ^ £l,i 

S{e,;x',y') = [ d'x ViT.(x',a;) / d'y V ^^^^^l^^^S^ Vi(2/,y'), (46b) 
Jl f^, Jl ^ - sl,i 

and, with the notation T'*™ (r, r' =t, J,, 5 = 0, 1, and m = —S, . . . ,<S) for the spin operator projecting the doublets 
T and r' onto their singlet [S — 0) or triplet {S = 1) combinations,^^ 



Aie.;x',y') = / d^x^^UA^^x) f d'y ^ "^■-^'^^'^-'^^^ (46c) 

Am{e.;x',y') = \vl:;, [ d^x V^,(x',a;) / d^y CAv.v') (m = 0,±l). (46d) 

^ Jl —f, Jl I ^v- eL,i 
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While the projector in Eqs. (46c) and (46d) ensures singlet and triplet symmetry for A and A^, respectively, 
the sums over the spin components a and a' involve all combinations of the eigenvectors and , i.e., contri- 
butions from both the singlet and the triplet anomalous densities. The product u£-|.(x)w£* (y), for instance, which 
Eq. (23d) associates with Xi(y, a^), contributes to A and to all three components of A^; should a global rotation of 
the quantization axis diagonalize both t{x' ,x) and t{y,y'), its contribution to the singlet potential A would vanish, 
so that the mismatch between the symmetries of the densities and of the KS potentials would be restricted to the z 
component of Xm ^m- In general, however, in the presence of magnetically inhomogeneous interfaces neither the 
z component nor the total spin on one side of the interface are conserved. 

In particular, if a singlet superconductor in region L is coupled to a normal metal in region R, the singlet anomalous 
density will contribute to the right-hand side of Eq. (46d). We expect, therefore, triplet Cooper pairs to be proximity- 
induced in the normal metal, along with singlet pairs due to the potential A. 

To conclude this section wc note that the coupling between anomalous densities with different symmetries is in 
line with spin conservation. The off-diagonal terms in our phenomenological Eq. (44) arise from magnetic degrees 
of freedom in the barrier, which interact with the spins of the electrons in L and R. Only the total spin = 
{Sl + St + SrY , which includes the contribution St from such degrees of freedom, must be conserved. While a 
tunneling matrix (44) diagonalizable by a uniform rotation of its spin variables is sufficient to conserve S\ and 5^, 
and while at least in special situations, this condition proves not necessary,**^ the right-hand side of Eq. (44) will in 
general fail to commute with S\ and hence allow <S ^ <S ± 1 Cooper-pair spin transitions. 

VI. CONCLUSIONS 

We have discussed a density-functional formalism describing superconductivity and magnetism, a generalization 
encompassing DFT for superconductors^'^^ and DFT for antiferromagnets^ that covers systems with coexisting order 
parameters. The new formalism being particularly practical in descriptions of two quantum liquids coupled by 
tunneling, wc have derived an exact expression relating the proximity field in region one to the eigenfunctions and 
eigenvalues of the unperturbed KS equations in the opposite region. 

The solution of the self-consistent cycle of equations presented in Section III B gives individual access to the singlet 
and the triplet Cooper pair densities, as well as to the staggered magnetization. This feature of the formalism was 
explored in Section IV B, which demonstrated that the proximity field induced by the singlet (triplet) density affects 
only the singlet (triplet) pairs and showed that DFT opens attractive perspectives for the study of the competition 
between order parameters in junctions such as a conventional or unconventional superconductor coupled to a non- 
collinear antiferromagnet. In Section V we considered the special case of a magnetic interface. We found that a 
singlet superconductor can in this case induce triplet pairing in a normal metal. We hope that our results stimulate 
further experimental scrutiny of macroscopically inhomogeneous quantum liquids and serve as the basis for an ah 
initio description of such systems. 
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